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PREFACE 


Research was conducted on the effect of elasticity in 
the construction joints on the auto-oscillation of spacecraft 
with gas-reactive direction systems; a mathematical model was 
obtained, on the basis of the method of mixed coordinates, of 
a generalized flexible spacecraft at one end of which was 
appended the directive action of a system of gas-reactive 
nozzles. Various structural forms were obtained functionally 
describing flexible spacecraft, as systems consisting of a 
solid central body with flexible structural elements joined 
to it. Studies of the auto-oscillatory processes were conducted 
on the basis of a method of point-by-point transformation. The 
work derives the equations of the correspondence function and 
the equations of velocity at the limiting cycle, taking into 
account the delays in the relay and command elements. These 
equations were studied using analog and digital computers. 
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EVALUATION OF THE EFFECT OF ELASTIC JOINTS 
ON THE AUTO-OSCILLATION OF SPACECRAFT 
IJITK GAS-REACTIVE DIRECTION SYSTEMS 


G. G. Sasin, 

USSR Academy of Sciences Institute of Space Research 

Both the scientific and experimental requirements of /5* 

contemporary artificial earth satellites (AES) have become 
more rigorous every year, as well as those with respect to 
the accuracy of stabilization. So, for example, the AES 
placed in geostationary orbits must be light and compact but 
at the same time must have a substantial lifetime. The greater 
and greater energy consumption of such AES has necessitated the 
use of large panels of solar cells (SC) . For the triaxial 
stabilization of a similar spacecraft at least one type has 
used systems of reaction motors of both hot and cold gases 
and will soon use electroreactive motors (ERM) . The wide 
frequency range of the spectrum of these engines will generate 
a variety of modes of oscillation. The effect of these high 
resonance frequencies on the characteristics of the roll 
regulator, the requirement of relatively sensitive reception 
in the 0-10 Hz band, the achievement of a high precision of 
orientation and stabilization (from several minutes down to 
seconds of arc) , can all lead to substantial difficulty in /6 

setting up such systems. 

Experience has shown that, with the growing requirements 
for accuracy of orientation and the high cost of mockups, 
there is an increasing value in many projected developments 
in creating a mathematical model of the spacecraft with 
heightened flexibility and complicated dynamics of construction. 


*Numbers in the margin refer to pagination in the foreign text. 
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since the performance of tests on the ground of such structures 
is rather expensive and often simply impossible. At the pres- 
ent time there already exists [1] an entire group of mathematical 
models describing flexible spacecraft, the most general outlines 
of which are that they are formed of an essentially solid body 
with elastic consoles attached to it, which in turn can be 
composed of a continuous medium [2] or of a discrete collection 
of adjoining solid bodies joined by means of massless elastic 
elements called finite elements [3]. Depending on whether the 
initial model of the panels is continuous or discrete, its 
displacement with respect to the parent body is characterized 
by either distributive or modal coordinates, while the orienta- 
tion of the basic body of the spacecraft in inertial space is 
always described by the three angles of its spatial position. 
Reference [4] proposed a method, combining the advantages of 
the discrete and distributive coordinates, called "the method 
of mixed (hybrid) coordinates." This method makes it possible 
to provide a mathematical description of the basic solid body 
with attached flexible panels in such a manner that the 
oscillation of the basic solid body is described with the 
usual differential equations and the oscillations of the 
flexible panels with equations in partial derivatives. 

In this work the method of mixed coordinates is used to /7 

obtain a mathematical model of a generalized flexible space- 
craft, to the basic body of which the directing action of a 
system of gas-reactive nozzles is added. The various struc- 
tural forms of the transmission function are derived, as well 
as the relationships among them. A method of pointwise trans- 
formation is used to study auto-oscillatory processes in the 
spacecraft, whose direction circuit is pictured as an aperiodic 
reverse signal (ARS) with a double time constant which includes 
a relay element with a band of insensitivity. The equations 
of correspondence which are obtained as well as the equation 
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of velocity in the limit cycle take into account the delays 
in the relay element and control console. The research made 
use of an analog computer of middling size (AVM-MN-18M) and 
a small digital machine (type MIR- 2) for which programs were 
developed which made it possible not only to study the transi- 
tional and established processes from the initial system of 
equations, to construct the Konig-Lamery diagrams, and to 
obtain the value of the velocity in the limit cycle, but also 
to maintain a double control on the correctness of the solution. 

Description of the Mathematical Model /8 


The dynamical model of the system is sketched in Figure 1. 
It consists of a solid, rigid body B and, attached to it, the 
flexible structures A, consisting in turn of the elementary 
bodies A^ . Assuming the elastic deformation to be negligibly 
small, let us first derive the equations of the superstructure. 

The Newton-Euler equations describing the motion of the 
element A-^ can be written as 

F l ~ mi ci pi 

f - H* 

In the inertial coordinate system the acceleration , Cl p: 
can be obtained as twice the differentiated sum of the dis- 
placement vectors X* C * & * 2* * U 1 , which determine the 
relation between the immovable point O' in the inertial space 
and Pi (see Figure 1) . 

Let the orthonormal system of vectors * it , refer to 

body B, and Q t , Q 3 ,Q, to body A. Then the transition 

from the system 5* to the system (where a = 1,2,3) can be 


( 1 ) 

( 2 ) 
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performed by means of the time-varying matrix of directional 
cosines, C. That is. 


■4 


...C/i Cfg C/3 


[41 

o 9 . 

- . 

£■>/ Cii C»s 

* 

l. 

s 


C31 C n C33 


w 


or i 

(°h cfg] (4) 

Geometrically, the coefficient is the cosine of the angle /9 

between the base vectors and the rotated, base vector 
f t &*& • We then write 

/ • _ _ 

QPi**- r ix+c*K+v+a t ) ( 5 ) 

1 7 T T 

The relation between the inertial base vectors <// * J 2 • and 
the base vectors t..i. J4 is determined by the matrix of 
directional cosines 8 , or 

ffl - 8{J] < 6 > 


Realizing that 


! %c- ‘<*/ dl 5 + 


we obtain from equations (5) and (1) 

pm OO 004 O O , * ! w ■» I - « 

where ( ° ) is the vectorial derivative in the system at rest 
with respect to the body B. 


The equation of rotational motion of element is 

fi= d/ d{ (d ! 'U) 1 ) •tf-i2 i +u> t xA t -u) t . ( 8 ) 
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Assuming a condition exists of small values of the 


orthogonal displacements 


A' • fit .A' of the flexible structure 


A relative to the body B, we obtain the inertial angular 
velocity of the element Aj_: 

U) ! - U>+J§‘ . 


(9) 


The vector lies within the base vectors l.e.j,. 
stituting equation (9) into equation (8), we obtain 


Sub- 


f ! m *js*+ *) • d* ub/li/S't/f'/jlt 


(10) 


Replacing u) with and with 

■Wjr. in equation (10) and ignoring 
obtain at last' 


we 


’ ■ ' O m 


(11) /lo 


The vectorial product of two vectors with the same coor- 
dinate bases is easily expressed in matrix form using a 
diagonally symmetric operator. For example, if. u) is any 
diagonally symmetric linear operator in the three-dimensional 
Euclidean vector space; , it is represented in the orthonormal 
base coordinates .1.1 by a diagonally symmetric matrix. . 
Thus , - - - - - 



0 U9, 

*> 3 O - «/y 

-u> 2 0 


( 12 ) 


and for each vector \£ within 'll , the relationship 

u)x c u)°C 


(13) 


is fulfilled. 
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Equations (7) and (11) can be rewritten in matrix form: 


f r -/n ! {8X+ e * u‘+2u)ft+u l ]+ T*u : ]+ u)-<2fc*frv*u‘j} ( (14) 


It is obvious that the matrix C which characterizes motions 
in the CM system depends on the displacements of the elements 
of the structure A, which are expressed in the matrix U 1 . 

But since, generally speaking, the matrix C also depends on 
other variables (the deformation of the other superstructure, 
external disturbances, directional moments, etc.) as well, the 
substitution for C should be written in the form 

C-^rtm : u : *e (16) 

I ' 1 i‘£ 

where e is any variation in the CM independent of the variation 
of the structure under consideration. 

Let us join the matrices in equations (14) and (15) . They 
have dimensions of 3x1 for each element A 1( . . . , A n in the 
simple matrix equation of the dimensionality 6n x 1. A total 
displacement motion (translational and rotational) of the 
structure is described by means of the matrix coordinates ^ 

(a matrix column of dimensions 6n x 1) . 

f = f uiu!ui/3}&&uteu;...p?] T (17) 

The inertia matrix y, which is zero except for the matrices 
of the mass and the moments oif inertia of the elements^ of the 
superstructure distributed along the main diagonal, has the form: 


/II 
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( 18 ) 




I. 


m ' f o 

rrf 

\ 

o r 

m 


Since the total matrix equation must contain elements 
which are characterized by two parameters ( W , and n l , I { ) , 
it is appropriate to introduce the matrix operators 

Z f0 =t [EQ£O...EQ] t (19) 

and 


z„w[om...oa T , (20 

where E and O are the unit and zero matrices, respectively, of 
dimensions 3x3. 


The matrix C in equation 

f » — 

M 


(16) now becomes 


llo *e 


( 21 ) 


We obtain the equation of motion of the structure, assuming 
the system is without friction or plastic deformation, from 
equations (14), (15) and (21): 

* 9J(‘E-?Sui -{e^P/ui * U>cd(e'f3j] + ( 22 ) 

where K is the matrix of the rigidities of the structure, 

A is the matrix of external forces and moments (of dimensions 
6n x 1) , and the new matrices ffj t $ ^ and i are defined 


9, s T fe cJ={e*'06*r...O] T 
9, s Z, t t*>*Cou>o...urT 
2 s Cl t OV...O] T 


(23) 


/12 


7 



(The tilda sumbol (~) above a matrix indicates dimensions of 
6n x 6 '. ) 


The extended matrix has elements of 3 x 3 dimensions and 
is zero off the diagonal, where the elements are the corre- 
sponding elements of diagonally symmetric matrices obtained 
in a manner similar to (12) and correspond to the matrices 
of dimensions 6n x 1. For example, obtained from equation 


(23) is 



and S2 t , Sl t , 9 t ,Q a 


and 


/-'S’, 


are similar. 


(24) 
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The Equation of Motion of the Spacecraft 


The equations of motion of the spacecraft are written with 
consideration of the fact that the flexible superstructures are 
contained in the structure of the spacecraft. Thus 

f-H 

■ * Mctt k fpxp dm 

We rewrite (25) taking (26) into account: 

where g 4 

We now rewrite (27) in matrix form: 


(25) 

(26) 

(27) 
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( 28 ) 


-fy 7 /* r * -T tt e r + 

■■■ -r t » - <>/« .. . . rimir. I T !/• 


1 - i'Bf’/tieo - ZeS/Sf. # 7 } +/2f/t(Z t 

7/2^ Tec- # T J “Sy' 


Linearization of the Equation of Motion 


With a wide variety of spacecraft which have flexible 
structures attached and are stabilized in 3 dimensions, it 
can be assumed that their angular velocity is small enough to 
ignore expressions of the second order and their derivatives. 

1 CO — Q ''W-- 

, & — £ 

6 — (£-81 

where a yre,B,e. T of the base vectors 4r»'4L4. connected to the 
body B. Under these conditions the linearized equation cor- 
responding to (22) is ' 

\-(e+R) 6] Tec B+A+/'t*> dX ( 29 ) 

We transform (29) according to the form of oscillation of the 
craft itself, 

• = <t>*l (30) 

* » 

and, rewriting (29) taking (30) into consideration, 


’-tW: 

-Tje>i?h ii n }B+rj(s-§li‘e- si a ]} - 
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Normalization of the proper vectors can be performed in 
the following manner. Let us assume that 

(where E is the unit matrix) . In that case becomes the 

matrix £ 9 , determining the proper values. These operators, 
in combination with equation (31), make it possible to obtain 
the equations of motion of a modular element in a form which 
in its homogeneous parts corresponds to the classical formula- 
tion of the equation of oscillation in the functions of 
normalized coordinates. 

In developments of the model of structures with discrete 
parameters the damping is usually ignored, but the diagonal 
damping matrix with the diagonal elements is included 

in equation (22) after transformation according to (30). 

In finite form the equation of motion of the flexible 
structure can be written in the form 

fa. ( 32 ) 

The linearized equation corresponding to equation (28) can be 
rewritten in view of the above premises, equation (30), and 
the fact that in small displacements of the structure the values 
of C are small enough to ignore expressions of the second order 
and their derivatives: 

T*IB TV <*>$( *■ TZ 2 +?«) (33) 


In order to choose the N modes of oscillation which will 
be appropriate and to represent with a sufficient degree of 
precision the dynamics of the reaction of the system which 
has been modeled by means of 6n modes of oscillation, certain 
criteria are necessary. For this criterion it is appropriate 


/15 
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to select the lowest modes or those modes in which the corre- 
sponding frequency will be located in the neighborhood of the 
assumed "forced frequency." 

Transmission Functions 


The equations of motion of the spacecraft which include 
the attached structures, (32) and (33) , after excluding the 
effects of external forces and displacements, and independent 
of internal displacements of the structures (that is 6 - 0 t 
<p r A - 0 ) , can be given in the following form: 


■ . f c <- 



(34) 

5 tPjfd iffy -a 8, 

(35) 

where 4 " 

Here and subsequently we regard the matrix of modal coordinates 
n of dimensions 6n x 1 as truncated to the f orm Q/ with dimen- 
sions N + 1, and the underline ( ) will be omitted. 

By using a Laplace transformation we exclude from 
eration the values of the displaced superstructure ri . 
equations (34) and (35) become 

consid- 

Then 

fysy-sfyeisi-s^Qis) 

(36) 

♦ ffyts) 

(37) 


Solving equation (37) for n(S) and substituting this into 
equation (36) , we obtain 

JT(S) •/£*/- S 4 A r m)d(S) ( 38 ) 

The matrix eis*6?r‘ is diagonal. In accordance with 
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the principle of superposition in the linearized system under 
consideration, we can rewrite the matrix of the transmission 
function of equation (38) in the form 




S'WAS+Cf 



(39) 


If the matrix of the transmission function is diagonal, the 
reaction of the spacecraft can be represented in the following 
system of equations: 


'rii *WQj * =■ Au6j+J ia § a 

Iidi+IitBg *I ia Q a +fk t /2 i m 

**T S 


A structural format corresponding to equation (40) is laid out 
in Figure 2. 


The scalar transmission function of the open contour 
G a (S) for the a axis (a = 1,2,3) is obtained from equation (39). 
When •Igj “Igy-Igt a J a g - 0 the transmission function becomes 


r /«_ SM • 


(41) 


The value in the round brackets is called the "normalized 
reduced inertia" for the a axis and the i th mode of oscilla- 
tion. It is designated /%J or simply R. 




The physical significance of R for the scalar values ^ 
or in a more general view for the (3 x 3 ) -dimensional matrix 
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T w T 

in equation (39) proceeds from a consideration 
of the limiting case. 

If the superstructure approaches rigidity, the transmis- 
sion function of equation (39) is reduced to the transmission 
function of a rigid spacecraft 

G($) = i/l A 'S* ( 43 ) 

If, on the other hand, the superstructure becomes elastic to 
the limit, we can assume that the transmission function will 
be based only on the transmission function of the solid body, 
while the superstructure will have essentially separated, 
since the reaction under consideration will be lacking on the 
part of the superstructure. Therefore, in that' case, the 
proper frequencies of the structure A in equation (39) will 
all tend toward zero, and the equation becomes 

! eis) -tfi-fc £njs‘]-*Tis) (44) 

That is, the 3x3 matrix should be the matrix of 

the moments of inertia of the spacecraft relative to its center 
of mass, without the moments of inertia of the structures. 

This value varies between zero and one, and the value from 
equation (42) can vary from zero !(for the limit of large 
structures which are very sensitive to deviations in &*) up 
to unity (for small structures or those not sensitive to 
changes in 8 * ) • Thus, 

In many cases it is more acceptable to write the trans- 
mission function describing the dynamics of the rigid space- 
craft in some other form than (41) . 


/18 
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Let us write equations (34) and (35) in scalar form, 
taking into account the fact that cross-connections between 
channels are lacking: 



(45) 


(46) 

where # - V& , 4* tf/i . 



From equations (45) and (46) we find 

JLrjt =-xi . 

Li Lt 


Designating 


we obtain 1 


/-L?Lt „ A . 

i WElT " ^ 

fa-Xi+fAtrii 
q] +2$6ift + Oi >, 7} m ~X 0 


The transmission function of this system will be 


r'/ol flffj '?U+Aih?X'S;'S'tf 
bJS)- TJS ) " pWVpWT 


The structural flow diagram corresponding to the transmission 
function of (49) and (50) is sketched in figure 3. To obtain 
equations (48) and (49) in dimensionless form, let us intro- 
duce the dimensionless time t° : 


p*,. — * - ms 

■ tFtfL? 

dc * r U i - dt = 6 'leu 
fFffL? 

de 2 * -rA=r cU‘- Sidi* 


Using (51) equation (47) becomes 
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(52) 


Al. dl ot df Ai U 


-x: 


Introducing the new variable 


At % ^ 


we obtain finally for the flexible element 




where 


$)* - ?/i fFtfEi 


Let us determine the expression for the coordinates of 
the solid body. Passing to dimensionless time t° in equation 
(45) , we obtain , 

d a 9*__ 4? {-tfiAt w . £ tr' d*!2i 


d a 8+ _ 4" At At "y* . sf AJ'JOllli 


Transforming this expression and using the relation (53), we 


obtain: 


4 F,' _ O . $ iTa, # 


Vie introduce the new variable 


* r & JL 


In final form the equation (55) becomes 




We will study the auto-oscillations of a nonrigid space - 
craft in the direction circuit of which we will include aper - 
iodic feedback. 


Description of the System 
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The system (see figure 4) includes a standard regulator 
with aperiodic feedback with a double time constant, a pair of 
gas-reactive nozzles, a position guage, and the dynamics of 
the nonrigid spacecraft. The regulator consists of a relay 
element, including aperiodic feedback with a double time 
constant (T», when ; yM-'S , r, when j frt) ’0 ) and an amplifi- 
cation coefficient -K 0 6. The zone of insensitivity of the re- 
lay element is -6-, the hysteresis in H. The gas-reactive 
nozzles have a thrust level F and a moment arm L. The dynamics / 21 
of the spacecraft include the oscillations of the spacecraft 
and the leastic weakly-damped oscillations of the attached 
structural elements. 

For the class of spacecraft under consideration — apparati 
with elastically attached structural elements — the initial 
synthesis and classification of the parameters of the circuits 
with aperiodic feedback is based on the premise that the 
spacecraft is basically governed by. the law of rigid bodies. 

The regulator should be designed in such a manner that the 
frequency of the proper oscillations of the superstucture will 
be above the frequency band of the channel of the direction 
system. It is however impossible to ignore the nonrigidities 
of the system: the damping of the proper modes of the elas- 

tic oscillations is very small, there are time delays in the 
circuit which are comparable to the periods of elastic modes 
of oscillation, the position gauge follows not only the low 
grquency of the motion of the solid part of the spacecraft 
but also the higher frequencies of the attached elements. 

Retruning to the previous section, we have obtained the 
equations of planar motion of a nonrigid spacecraft in dimen- 
sionless form, namely: 
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( 58 ) 


where 


& ♦*<*>/ 

(K.S)f~z+(r:,i)5 t 



* X npu LL>S 

0 npu -£S<U<£' 

- ^ npu U< ti 

' . a i 

- X npu u*-f 

0 npu -£<U<€$ 

* npu U 


for U 


for U<0 

. : ' . .1 


f 



tom npu 


f U~SuU>0 

\U°-6j 4 Lt*0 

fu r ?£ u u<o 
[u-efuu>o , 


(54) 

(59) 


(60) /22 


(61) 


Equations ' (59) , (60), (61) describe the action of the regulator. 

The system has delays in the relay element and in the 
command element 1 tj. A flow diagram of the stabilization circuit 
of a nonrigid spacecraft is presented in figure 4. At the in- 
put of the relay element (1) the directing signal U(t) arrives; 
at the output of the relay element the signal f (t) enters the 
command element (2) . The command element generates the stabil- 
izing moment T(t) which acts on the spacecraft in the desired 
manner. This same signal f(t) includes feedback (4). The 
feedback signal Z (t) and the output signal of the angle gauge 
(5) enter the summing device (6) , in which the directing signal 
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U(t) is formed. 


We assume that disturbing moments do not act on the sys- 
tem, that the angle gauge is ideal, and that the command and 
relay elements used are ideal circuit elements with pure time 
delays. The nonrigidity of the system is incorporated into the 
model as an externally attached elastic structure (7) . 

The real motion of the dynamical system corresponds to 
the motion of a representative point in phase space with the / 23 

coordinates b . d,. i . Let us project the trajectory of the 
motion of the representative point onto the phase plane B .6 , 
passing from the consideration of spatial motion to planar mo- 
tion. 


We find the solutions of equations (58), (54), (59) in the 

mth section of the trajectory, where % * Q : 

a) for the successively integrated equation (58)., we obtain: 


fff- •k-iAifa 


(62) 


b) for equation (54) the solution has the form: 

* ^ u sln W , (63) 

hr cos w- 

where vi* i irw and )^ #1 are the initial conditions. 

c) the solution to the third equation has the form: 


iay— K.8(i-e' t/a ) 


(64) 


We begin consideration of the motion at the moment t° = 0, 

when 
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I ®!i^0 


^ I I Vo 0 ' 


The command element is switched on at the moment *7 l b . At 
this moment (see figure 4) 

o (65) 

Ate, 

Beginning from the moment , the representative point 

moves along a phase trajectory intersecting the axis , §-Q • 

At the moment in time , when 

ie-x (66) 

the signal for switching off the command element goes from the 
ouput of the relay element to the electropneumatic valve (EPV) 
of the command element. 


The equations of motion remain as in (62) and (64), but 
taking into account the delays they assume the form 


2(W+t'J,h-Kj u-e 

Bt* -Xtibt.y+FAi/ii Ai q t , 

r i’l 


Taking into account the fact that at the moment of each 

subsequent switching off of the command element the elastic 

oscillations of the flexible structures diminish, that is, that 
.. ' A . n 

ai*- a* - u . the latter two equations of system (67) can be 
rewritten, taking (65) and (66) into consideration, where 

sm-tv-s,’. ife-w-% ■ we have 


1C 
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Switching off the command element takes place at the moment 
in time • At this moment 


4- 4r£<}t 


(69) 


Taking (69) into account and assuming that; , 
(that is, that the delays of the switching on and switching off 
are the same) the system (68) can be rewritten 


where 


X d x * 

Qr-W+TWt+B- 


(70) 


/25 


ASc-h 

Noting that 


^S 1 and 


S in the system of equations 


lv 


(70) and solving this system for S and S', we obtain 

• ; fi ■ ' 

ri r -K>$u-e 


(71) 


where 


%— Zn-e** t ?caritG+ ] 

< ni - •%/$'£***< sin till 


(72) 


Beginning with the moment ** ■ <7 ♦ tj the representative 
point intersects in its trajectory the $ axis, and at the 
moment in time i* - tg * if , the signal to switch off the 
command element goes from the relay element to the EPV of the 
command element, which ceases to operate at the moment t° = 



Then for the segment where the directing moment disappears 
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-A 

( X -0 ) * the equation for the velocity of motion of the system 
becomes 

(73) 

A, 4 

where file “ 1 i from the second equation of the system (72) . 


Because the system under consideration has no unsymmetrical 
limit cycle, since in the idealized version presented here no 
disturbing moments act on it, the function of resenblance can 
be examined in the sections S, S', and S". In the final varia- 
tion, taking equation (73) into account, the function of corres- 
pondence can be written 

f 9/2r(; 3 -2iu;+t?)+Ko6(i-e~ t,/n )-H -t Ai ^ 


(£♦> 


„ v/pvf+pfi e,-t1\-KoSa-e <;/r V# * 'iMj 




ti 


fa?»*r4*e OSM- MJ.0 ; 


where 

are the conditions of attenuation. 


(74) 
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Determination' of the Velocity in the' Limit Cycle 


From the first equation of the system (70) we have 

_ s*s" 

1 ‘ = J __ J 


Substituting (75) into the second equation of system (70) and 
taking account of the fact that $ ~ § 0 • S" - - ' we 




(76) 


obtain 

2/i n"^rrr~i .( 

Tf ,J J 2 to 

From equation (76) we find the velocity at the limit cycle q { 
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( 77 ) 


for which S * = S ** 3 0 . 

'< K.ht-e^i-2i e, ♦ ^ ;*-£ m, 


or 


whence 




e -'x¥j s /_ H/K.8-2@<Sf<.h% t^KoS-IAtq/tC 8 , 


Uo ■ * ■&(* 7&-l&r+^--£F-}' 


( 78 ) 


( 79 ) 


Results of the Modeling 


The analysis described above was conducted for the stabil- 
ization of a spacecraft along a single axis. Conducting an ana- 
lytical study of the stailization of the spatial motions of a 
spacecraft represented a task of great complexity, so at the 
design stage the real direction devices were designed to be 
replaced with modeling, using an electronic computer. In addi- 
tion, at the planning stage the interactions determined by the 
values of I a , I a > # Ijt • 1st • (see figure 2) could frequent- 

ly be ignored, due to the symmetries of the system and the small 
angular displacements in the area of stabilization. 

The modeling problem involved the study of the dynamical 
behavior of the system for one total moment of inertia I of the 
system but different variants of the value of the moment of iner- 
tia I of the solid body and I of the attached elastic structural 

K n 

elements (panels of solar cells) ; for different values of the 
frequency of proper oscillation of the solar cell panels; with 
and without time delays. We also verified the possibility of 
performing similar studies both on an analog and on a small 
digital computer (of the Mir- 2 type) , for which programs were 
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developed. For the solving of the equations on the analog com- 
puter (type MN-18M) the structural model of the transmission 
function was obtained. 


Equations (58), (54), (59), and the corresponding flow dia- 
gram in figure 4 were modeled in the analog computer in the form 
of the chart in figure 5. The delay was modeled according to the 
Pade arrangement (using amplifiers 10, 11, 12, 13). The elas- 
ticity, including the fundamental frequency of elastic oscilla- 
tions, was modeled either by the amplifiers 6, 7, 8, 9 (which cor- 
responds to the structural arrangement in figure 3a) or by 
amplifiers 14, 15,16,17 (which corresponds to figure 3b). The 
selection of the elements R^, R^ , of the aperiodic feed- 
back with a double time constant was made in the following 
manner. When the circuit is turned on, the signal Z (with 
zero initial conditions) at the output is 




(80) 


At the switching off with Z 0 initial conditions the signal Z 
becomes 

zut=z 0 e~ t/ * ,c * 


(81) 


So that 
or 


fi • • RiCf » To = RiRzCg/(R t +Rg) 
R t .** Tt/Ct , Rg * % 


In the determination of the velocity 6a in the limit 
cycle equation (77) is used, transformed to the form 

X. t? 


ti 'll, 


(82) 


where 




( 83 ) 


and 
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The given equations were solved in the digital computer by 
means of the program in figure 8. In order to solve the same 
equation on the analog computer, a Laplace transformation was 
performed, obtaining 



(84) 


The model with the transmission function (84) is shown in figure 

6 . 


In order to construct the Konigs-Lamery parametric equa- 
tion, the correspondance fuction (74) is rewritten in the form 


+ K*SU~e~ it/tQ )-H - i A:ai 

^ as X, (i\ u -$)-K a 8U-e ' u/ta bH + £ Mi 


(85) /29 


This equation was solved on the digital computer by means of the 


program in figure 9. 

sn); 


The printout produced the values of t-^, 


where 




( 86 ) 


The structural model of the system (85) for the use of the 
analog computer was obtained in a manner analogous to (84) 
and is diagrammed in figure 7. 


The parameters of the system of equations and the dynamical 

characteristics of the model under study are produced for various 

values of I and I in Table 1 and presented in figure 9. The 
K n 

same system was studied with weak damping of the proper modes of 
elastic oscillation with the following parameters: 
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sec 


K 0 S = 162.9 arc min.; H = 0.9 arc min.; x 0 = 2.1 sec; =39.6 
T = 0.682 N ; I = 1115 Nm sec 2 ; Xo = 2.1 arc min/sec 2 

A =6.39 ; 44 = 803.9 8 Nmsec 2 ; 4 =0.1 Hz 

A = 0.853 ; A4 = 107.3 Nm sec 2 ; A = 0.18 Hz 

A = 0.457 ; AA = 57.56 Nm sec 2 ; A‘= 0.224 Hz 

?/•?*" = 3.5 X 10" 3 , ‘ 

It was shown in the modeling process that it is possible 
to obtain solutions to equations (84) and (85) on both an ana- 
log and a digital computer. 

The solution to equation' (82) for the 20th variant of Table 
1 is given in figure 10. 

The numerical solution to the system (85) is given both 
selectively in Table 2 for the variants of the 1st, 5th, 7th, / 30 
9th, and 15th modes of the Konigs-Lamery diagram including 
delays ( t s = 0.1 sec for "1" , "5" , "9" , "15" and tj = 0.2 sec 
for "1") and in figure 11 without them. 

Figure 12 shows . the dependence of the velocity in the 
limit cycle , As on the inertial coefficient characterizing 
the flexible superstructure A with time delay. 

We should note that upon an increase in A the velocity 6U 
decreases, and in the presence of time delays it declines par- 
ticularly rapidly as A > 10. As A ■* 0 the effect of a time 

_ ^ 

delay on the change in Qe is insignificant. As is clear from 
an examination of figure 10 and Table 2, the change in value 
of the frequencies of the proper oscillations of the super- 
structure, 0* with a single value of A, are practically un- 
• • 

affected by - Q a , especially for A < 10. By means of the 
modeling on [illegible], shown in figure 5, confirmed the 
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complete identity of the structures modeling the dynamics 
of a spacecraft with attached flexible structures in Figure 3. 
The research conducted on the analog computer showed that the 
largest contribution to the dynamical processes of the system 
belonged to the primary (first) harmonic of the elastic oscilla- 

r. 

tions of the structure and that at the initial design stages 
of the system, higher harmonics can be neglected. 
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[*Commas in tabulated material are equivalent to decimal points. 



TABLE 2. 


I”10 

"TAJJMU'i 


I 

S 

S? 31 

3H 

.5000,-1 -.1100, 1 

.2150b 1 -.1579.-1 

.1051. 10 

.3020.-1 -.1797, 0 

.1236. 1 -.1600.-1 

.105fb 10 

.3040.-1 

.7323. 0 

.3316,0 -,1622.-1 

.106^, 10 

,306Cb-l 

.1635. 1 

-.5649. 0 r-. 164 3.-1 

,1072. 10 

.3000.-1 

.'2531. 1 

-.1453. 1 -.1665,-1 

.1079. 10 

.3100.-1 

.3418. 1 

-.2333, 1 -.1686.-1 

.1086, 10 

.3120b-l 

.4297. 1 

-.3205. 1 -4708,-1 

.1093, 10 

.3140b”l 

.5168, 1 

-.4069. 1 -,1730.-1 

.1100. 10 

.3160.-1 

.6031. 1 

-.4925. 1 -.1752,-1 

.1107, 10 

.3180b-l 

.6006. 1 

-.5773. 1 -.1774.-1 

.1114, 10 

.3200b“l .7735. 1 

1-20 

TAU)IHUA*1 

-.6615. 1 -.1797.-1 

*1121. 1 

T 

S 

S2 01 

3H 

.3000.-1 -.2082. 1 

.3132. 1 -.1540,-1 

.1012, 10 

.3O20b-l -.1174. 1 

.2231. 1 -.1560.-1 

.1018,10 

.3040b“l -.2754, 0 

.1339b 1 -.1580b-i 

.1024, 10 

.3060b A l 

.6148b 0 

.4561. 0 -.1601.-1 

.1031. 10 

.3080.-1 

.1496, 1 

-.4187, 0 -.1622.-1 

.1037. 10 


.3100b-l .237Gb i -.1285. 1 -.1642.-1 .1004. io 

.3120.-1 ,3235b 1 -.2143. 1 -.1663.-1 .1050b 10 

.3140.-1 .4093.1 -.2994,1 -.1604,-1 .1057,10 

,3160,-1 .4942,1 -.3836. 1 -.1706,-1 .1063,10 

.31«V*1 .5704, 1 -.4671. 1 -.1727,-1 ,1069b 1° 

.3200^-1 . .6610b 1 -.5490, 1 -,1740b-i .1076, 1 

I“50 

,3900b-l ..2717,1 -.1352.1 -.2653.rl .1359,10 

.3920b-l .3210b 1 -.1830.1 -.2600.-1 ,1366b 10 

s .3940.-1 .369^ 1 -,2320b 1 -.270V1 .1373b 10 

.3960.-1* .4184, 1 -,2790b 1 -.2735.-1 .1380, 10 

.3980.-1 .4666, 1 -.3273. 1 -.2763,-1 ,1387, 10 

,4000b-il .5143, 1 -.3743. 1 -.2791.-1 * .1394. 10 

.4020,-1 ,5617. 1 -.4210, 1 -.2819,-1 .1401. 10 

.4040.-1 ,6(W, 1 -.4673. 1 -.2047.-1 .1400. 10 

.406 Ob -1 .6553b 1 -.5132. 1 -.2875.-1 .1415. 10 

,40«0b-l .7016,, 1 -.5588. 1 -,2903.-1 .1422. 10 

.4100.-1 .7475,, 1 -.6040. 1 -.2932.-1 .1429. 1 
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TABLE 2 (continued) 


I**40 


. 3900.-1 

-.2455^.1 -.1090,, 1 

-.2610ft-l 

,1332ft 10 

.3920^-1 

.2945 tt 1 -,1573 b 1 

-,2644 u -l 

.1339. 10 

.3940,-1 

.3431. 1 -.2052b 1 

-.2671 b -1 

.1345. 10 

. 3960^-1 

,3914b 1 -.2528b 1 

-.2698.-1 

.1352ft 10 

.3980^-1 

.4392o 1 -.2999 tt 1 

-.2725ft-l 

.1359b 10 

iAOO’Vl 

.4 (367b 1 -. 3467 b 1 

-.2753.-1 

. .13651s 10 

.4020.-1 

.5338ft 1 -,3931 b 1 

-.2780ft-l 

.1372. 10 

.4040^-1 

.5805b 1 -.4391b 1 

-,2808.-1 

•.1379, 10 

.4060^-1 

; . .6269b 1 -.404B n 1 

-.2035 s -l 

.1385, 10 

.4080.-1 

.6729ft 1 -,5301ft 1 

-.2863„-l 

.1392, 10 

.410Qb-i : 

.71G5n 1 -,5750* 1 

-.2891 D -1 

.1399. 1 

1*50 




: .4 30^-1“ 

.7270b 0 .7779ft 0 

-.3253.-1 

.1512„ 10 

\ .4310 o — i 

,9133*, 0 .5951b 0 

-.32o0ft-l 

,1516ft 10 

.43200-1 : 

.1098ft 1 .4131» 0 

-.3204 o -i 

.1519. 10 

•4330.-1 

.1283b 1 .2318ft 0 

-.3299.-1 

.1523^ 10 

.4340.-1 

J . 1467 0 1 .512<V1 

-,3314 b -1 

.1526. 10 

.4350.-1 

.1651 R i 1 -,1285» 0 

-. 33?9„-l 

.1530b 10 

.4360,0-1 

.1833ft 1 -.3076 b 0 

-.3345ft-l 

.1533. 10 

.43700-1 

.2015ft 1 -,4859 b 0 

-. 3360„-l 

.1537„ 10 

.4380.-1 

.2196„ 1 -.6636,, 0 

-. 3375.-1 

.1540,, 10 

,43900-1 

•2377ft. 1 -,8405 b 0 

-.3391.-1 

.1544. 10 

.4401)0-1 

.2556b 1 -,1016ft 1 

-. 3406 b -1 

.1547, 1 

1*60 




. ,4300.-1 

*, 6493n 0 ,0556ft 0 

-.3229„-l 

.1496. 10 

.4310,, -1 

.0352ft 0 ,6732 b 0 

-.3244.-1 

.1499. 10 

: .4320.-1 

.1020ft 1 .4915» 0 

-,3259„-l 

.1502. 10 

.43300-1 

.1204ft 1 .3106 b 0 

-,3274b-1 

.1506, 10 

.4340.-1 

.1380b l . 1304b 0 

-.3290.-1 

.1509® 10 

,4350.-1 

•1571b 1 -,4897ft-l 

-. 3305ft-l 

.1513. 10 

.4360„-l 

,1753b 1 -.2276 D 0 

-.33 20|j- 1 

.1516, 10 

.4370.-1 

. .1935ft 1 -.4056ft 0 

-,3335b-1 

.1520. 10 

.4300„-l 

.2115ft 1 -.5829 n 0 

-.3350.-1 

.1523. 10 

' ,4390b-1 

.2295b 1 -.7594n 0 

-,3365b-1 

.1526*, 10 

,4400b -1 

.2475b 1 -.9353ft 0 

-.3381.-1 

.153CV, 1 


nr v Slv * 
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TABLE 2 (continued) 


l»7o 


.4400.-1 -.1351. 1 

2091. 1 

-.3388,-1 

.1539. 14 

.4420.-1 -.1012. 1 

.2559. i'. 

-.3419.-1 

.1546, 16 

.4440^,-1 -.6755. 0 

.2229. 1 

-.3450,-1 

.1553, 16 

,4460.-1 -.3414. 0 

.1902. 1 

-.3481.-1 

.1560^ 16 

,4480.-1 -.1009.-1 

; .1570. 1 

-,3512.-1 

*1567. 16 

.4500.-1 ,3185. 0 

*1256. 1 

S 3544.-1 

.1574, 14 


.4520.-1 *6446*0 .9373.0 -.3575.-1 .1581,10 

.4540b-l .9680. 0 *6209. 0 -.3607^-1 *1588. 19 


' .4560.-1 

' *1288. 1 

,3070. 0 -.3639.-1 

.1595. 19 

*4580.-1 

*1607,1 

’-*4342.-2 -*3671.-i 

.1602, 16 

.4600.-1 

1-86 

.1923,1 

-.3132. 0: -*3703,-1 

1 ■ 

.1609, 1 

*4400b-l 

-.1388. 1 

.2928. 1: 3370,-4. 

.1527; 10 

.4420.-1 

-,1049. 1 

’.2596. i -‘,3401.-1 

.1534, 10 

.4440.-1 

-.7130^ 0 

.2267, 1 -*.3432.-1 

.1541, 10 

.4460.-1 

“43792. 6 

,1940^ 1 ♦•*3462.-1 

.1548, 10 

.4480,-1 

-,.4827.-1 

,1616. 1 -> *93.-1 

.1554, 10 

.4500.-1 

>2800. 6 

.1294, 1 -.3525,-1 

.1561. 10 

.4520^-1 

,6057. 0 

,9762. 0 -.3556.-1 

.1568, 10 

.4540.-1 

.9288. 0 

,6601. 0 3587]j,-l 

.1575, 10 

.4560.-1 

.1249. 1 

,3465b 0 -.3619.-1 

.1582, 10 

.4580.-1 

• .1567. 1 

.3558,-1 -.3651.-1 

,1589, 10 

.4600,-1 

1-96 

,1882. 1 

-.2729. 0 -.3683,-1 

.1596, i 

,.4500.-1 

-.1653, 1 

,322^1 -.3542,-1 

.1574. 10 

.4520,-1 

-.1336. 1 

.2918, X -.3573,-1 

,1581. 40 

,4540^-1 

-i.1021. 1 

.2610, 1 -.3605,-1 

.1588, 10 

.4560,-1 

-.7092, 0 

.2305, 1 -.36.37,-1 

.,1595, 10 

.4580,-1 

-.3996. 0 

.2002, 1 -.3669.-1 

ji602, 10 

,1609, 10 

.4600.-1 

-.9254,-1 

.1702, 1 r. 3701,-1 

.4620^-1 

,.21 20. ,0 

.1404. i. -,3733.-1 

.1 61- c u 10 

,4640,-1 

.5142. 0 

*1109. l! -,3766.^1 

.1623, 10 

,4660,-1 

. 8140* 0 

.0169, 0 -.3798,-1 , 

,1630, 10 

.4680,-1 

.1111. 1 

.5264. 0 -.3831.-1 

.1637, 10 

; .4700,-1 

.1406. 1 

.2383, 0 -.3864.-1 

.1643, 1 


I 

i 

i 

! 
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1-100 
.4500,-1 
.4520.-1 
.4540^-1 
.4560.-1 
; 4580 B -1 
.4600,-1 
.*620*1-1 
.4640,-1 
.4660,-1 
. 4660 * 1-1 
.47(JQ,-1 
l4l0 . 
.4600,-1 
.4620*,-! 
.4640,-1 
.4660.-1 
.4680,-1 
.4700,-1 
.4720,-1 
,4740*,-l 
,476<W 
.<*700*,-l 
.4800^-1 
1-126 
.4600*1-1 
. 4620 , 0-1 
,464C\j-l 
.4660^-1 
.4600,1-1 
.4700^-1 
.4720*1-1 
.4740,-1 
,476Ct-l 
.4780,-1 
.4'8U0,-1 


TABLE 2 (continued) 


-.1673b 1 

,3248, 1 -.3528,-1 

-.1564. 16 

-.1356,, 1 

,2938b 1 -.3560,-1 

.1571, 10 

-.1041*, 1 

.2630, 1 -.3591,-1 

.1576, 10 

-,7296b 0 

.2325, 1 -.3623,-1 

.1585, 10 

-.42U2*i 0 

.2023, 1 -.3654,-1 

,1592. 10 

-.1133b 0 

,1723, 1 -.3606,-1 

,1599. 10 

.1911* 0 

,1425b 1 -.371^-1 

,1606, 10 

.4931 b 0 

.1130* 1 -.3751,-1 

.1613, 10 

,7927b 0 

.0382, 0 -.7/8V1 

.1620. 10 

,1090b 1 

.5479, 0 -.3815.-1- 

,1626. 10 

,1364b 1 

.2600,, 0 -.384$,-i 

.1633, 1 

-,1291b 1 

.2901. 1 -.3708,-1 

.1612. 10 

—.9925, 0 

.2609, 1 -.3740,-1 

.1619. 10 

-.6955b 0 

.2319, 1 -.3773,-1 ' 

.1626. 10 

1.4009,, 0 

.2031, 1 -.3805,-1 

.1633, 10 

-,1087b 0 

.1746, 1 -.3838,-1 

.1640, 10 

.1011, 0 

.1463* 1 -.3871,-1 

.1647, 10 

.4667*1 0 

.1183,1 -.3904,-} 

.1654. 10 

,7540b 0 

.9049,0 -.3937,-1 

,1661. 10 

.1037*, 1 

,6209, 0 -.3971,-1 

.1668. 10 

.1317 B 1 

,3551, 0 -.4004,-1 

.1675. 10 

,1596b 1 

.8353,-1 -.4038,-1 

.1682, 1 

-,1306b l 

,2916. 1 -.3694,-1 

.1603, 10 

loft/, 1 

,2624, 1 -.3727,-1 

.1610, 10 

-,7106*. 0 

,2334, 1 -,3759„-l 

.1617. 10 

-.4161. 0 

,2047, 1 -,3791.-1 

.1624. 10 

-,1240b 0 

,1762, 1 -.3824,-1 

.1631. 10 

,1656 r 0 

,1479, 1 -.3857,-1 

.163$, 10 

.4531. 0 

.1198, 1 -.3889,-1 

.1645, 10 

,7302. 0 

,9207, 0 -.3922,-1 

.1652. 10 

,1021. 1 

,6448. 0 -,3955*-l 

.1659, 10 

1301, 1 

.3711, 0 -.3909.-1 

.1666. 10 

.1500. 1 

.9963,-1 -.4022.-1 

.1672. 1 
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TABLE .2 (continued) 


1**30 

.0700.-1 

,i»'/20b-l 

W 

B-l 
.0760.-1 
.0800,-1 
.0020,-1 
.0040.-1 

,006 Oo-l 
.0000.-1 
.09OCb-l 
WOO 
,. 0700.-1 
.0720.-1 
.0700^-1 
,0760.-1 


- 46935. 0 
-. 3695 . 6 
-. 8776.-1 : 
. 1917 . 0 
. 0669 . 0 , 

. 7000, 0 
.1016. 11 
.1207. l.;H 
'•.,1556.. 1 i" 
, .1022. 1,1 
.2007. i i 


.2290. i 
,2021. 1 
.1706. 1 
.1070., i 

4 l 2 oo. 1 

,9359. 0 
.6700. 0 
.0063. 0 
.1006. 0 
’,1100.0 
ttO; 


■ •‘,3850.-1 ;v 
r, 3891.-1 ; 
-,3920.71 i 
■->3957.-1 . 
-,3991.-4 
-.0020.-1 ; : 
-.0058,-1 •: 

. >.0092.71 t 
-,0125.-1! 
-.0159.-1 
-.0190.-1. 


,1601. 10 
,1608. 10 
,1655. 10 
1,1662.: 10 
: .1669. 16 
. .1676. 10 
4 16 83. 10 
, .1690. 10 
,1697. 10 
,1700.. 10 
.1711. 1 . 


• .14600^,-1 
,0820^1 
••0800b"! 
.0860.-1 
•0880,-1 
,0900b-* 

,.08pct-t 
V .0820^1 
.•:(-:-> : 080Ctffi 
M->08iGtrl 
•oeaot-i 
V:.0900t-i 
•0920W-1 
; ,0900.-1 
; 1096CV-1 
.09 80.-1 

lV.JPD0.-rl: 


. r.6623. 0 Si 

.2307. 1 -.3808.-1 

4635. 10 . 

7,3780. 0 

.2030. t 7.3881.-1 . 

.1601, 10 

7,9678.-1 !; 

,1755. 1 -.3913.-1 

,1608. 10 

•*1826. 0,1 

•1083. 1 -.3907.-1 

.1655, 10 

;,0598. 0 

>1213.1 7,3980.-1 

,1662. 10 

;,7307. 0 ; 

•9052. Or. 4013.71 

.1669. 40. 

I.lflOT. 4 ) 

>6793. 0 -.4007.-1 

4676. 10. 

;:i;,l506i:iv'i 

;;0157. Q • -.0080,-1 

: 4683. 10 

.1501. 0. -,0110.-1 

>1690. 10 

■••!:.1013W^H 

—,1052. 0 — ,0108,, -1 

4697. 10 

;-4;i)7:7 i.;!'::,. 
•• *- -}.■ ( 

7,^6. Q -,0182.71 

.1700. 1; 

7,9292. O ' 

.2609. i -.0033,-1 : • 

j ,1680b i0 

7,6633. 0 

,2350. * -40067.-4 

.1667. 10 


: ,2093b 4 r. 4100.-1 

.1690. 10 

> -*1379. o 

^1038b’4' : ,;-,4130 mT l- 

>1701. 10 

..•:.Vl2i6i-::bV 

•1586b ilr, 0168.-14' 

■ ,1708, 10 

3791. :0' :i > 

4 1335b 4 V-,0203.-i\ 

.>1715. 10 

[,•6306. 0 ;, 

;;j087b'i : -,0237.-4. 

v 4722. 16.' 

v'i 8880^ :(>;•,. 

:>.£0?7?bril ;: .'’ 

i 4729, 10 


,5960b 0 -vo 306^-4 

‘M736^ lb 

Air.1389.-jv! 

•rV3539AO:'7.034i:^J; : 

. ■1703, JO 


. 4113^; 0 >,4376.-1; 

4750. i‘:: 
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TABLE- 2 (continued) 


M tfi 



.1676* l4 

,40Q0*-1 

-.9321,, 0 

.2612* 1 -.4026*-l 

.4820.-1 

-,to j62* 0 

.2353* 1 -.4060*-! 

.1683* 16 

.4040.-1 

-. 4024 * o 

,2096* 1 -.4093.-1 

.1690. 16 

» 4060.-1 

-.1408. 0 

.1041* 1 -.4127.-1 

- .1696* 16 

•WtKV-1 

*1 187* 0 

.1509*1 -.4l6l„-l 

.1703* 16 

.4900.-1 

.mi* o 

.1330* 1 -.4193,-1 

.1710* 10 

.4920.-1 

.6316,, 0 

.1090* 1 -.4230.-1 

.1717. 16 

,4940.-1 

.0U5(V, 0 

.0439* 0 -.4264.-1 

.1724. 16 


.11 36* 1 

.5995* 0 -.4299.-1 

.1731* 16 

.A9UCL-1 

,1383* 1 

.3570* 0 -,4 333.-1 

.1738* 16 

.5000.-1 

.1633* 1 

.1165* 0 -.4368.-1 

.1745* 

1-176 

.4900.-1 

-.4805* 0 

;2195* i -.4193.-1 

.1711. 16 

.4920.-1 

0 

.1950* 1 — ,4227.-1 

.1718. 10 

.4940.-1 

.21 3 VI 

.1707* 1 -.4262.-1 

.1725. 16 

.4960^-1 

. 2693* 0 

.1466* 1 -.4297.-1 1 

.1732. 16 

.4980.-1 

.5152* 0 

.1227* 1 -.4331.-1 

.1739, 16 

.®00*-l 

,7593* 0 

.9906* 0 -.4366,,-t 

.1746* 16 

• 5020U-1 

.1(X)1 0 1 

.7503* 0 -.4401.-1 

.1753*, 16 

.5040.-1 

. 124.1 * 1 

.5223* 0 -.4436.-1 

.1760. 16 

.5060.-1 

.1600* 1 

.2909* 0 -.4472.-1 

.1767. 16 

.3080*-! 

.1716* 1 

,6145*-1 -.4507.-1 

,1774., 10 

.3100.-1 

. 1951* 1 

-.1662. 0 -.4542.-1 

. 1781. 1 

1*4 00 
.4900.-1 

**.4013* 0 

.2196. 1 -.4188.-1 

.1708*16 

.4920.-1 

-*.2293* 0 

.1951* 1 r-.4223.-l 

.1715. 16 

. 49 40.-1 

J2 037.-1 

.1700* 1 -.4257.-1 

.1722. 16 

.4960.-1 

.2603* 0 

f 1467. 1 -.4292.-1 

. 1729* 16 

,4980fc-l 

.514?,, 0 

. 1228. t -.4326.-1 

.1736* 16 

.5080.-1 

.7303* 0 

.9916. 0 -.4361.-1 

.174 3. 16 

.3020.-1 

. 1000* 1 

.7565.0 -.4396.-1 

.1750* 16 

.10 40.-1 

.1240* 1 

,523 3* 0 -.4431.-1 

.1757. 16 

.a)60^-i 

.1479* 1 

,2919* 0 -.4466,-1 

.1764. 16 

.3000^,-4 

.1713*1 

.6249.-1 -.4502.-1 

.1771. 16 

..100.-1 

.1930* 1 

-.1631. 0 -.457/*-! 

.1778. 1 



TABLE 2 (continued) 

'• " . ‘ "■ ' '' 

• t ; .4900*?* : -i-. 9060b 0 . .2621.1 -.4217.-1 .1721. 10 • 

.... ^.*920^ . r.6557. 0 ] . 2377. 1 7.4252.-1 , .1728. 10 
V;: .'940.-1 ^-.4075. 0 ' .2130. 1 *%42 86.-1. .1735. 10 

’• •& 60.-1 -.1613. 0 .1097. 1 -.4321.-1 .1742. 10 

* .. . iBeiVri .1660. 1 -.4356.-1 *1749.10 

.5000.-1 .3251. 0 .1424,, 1 -.4391.-1 .1756. 10 

-*■ • 10 20.-1 V l , 565 5. 0 • ;l 191. 1 -.4426.-1 , 1763. 10 

•” , 50 4Qt -1 ,8040. 0 .9599. 0 -.4462.-1 .1774 10 

.5064*1 .1044 1: 1 7303. 0 -.4497.-1 .1777. 10 

/ .50 00b-l .1275. V .5025. 0 -.4534-1 .1704. 10 

.51 OQb-1 . 1508. 1 .2766. 0 -.4569.-1 .1791. 1 1 

aeco.-l 19405. 0 .2655. 1 -.3879.-1- ' ,1582. 10 ’ 

: .10 20^-1 r*6 904. 0 , .2412. 1 -.3910.^-1 .1589. 10 

• Wititri -.4424. 0 .2171. 1 -.3942.-1 ,1595. 10 

. >9 60** 1 -.1963. 0 .1932. 1 -.3974.-1 .1602. 10 > 

; .«G0rf-l - .4776.-1 » 1695. 1 ;-i4 006.-1 *; .1608. JO 

.IDOQj-l .2099. 0 .1464 1 -.4039.-1 . 1615. 10 '■ 

v . 5024*4 ;5301. 0 I. ' , 1226. 1 "-.4071.-1 , 1621. 10 

C .5044-1 ,7605. 0 . 9954. 0 -.4104.-1 .1628. 10 

^ .50 60.-1 .1005. 1 .7659. 0 -.4136.-1 . .1634. 10 

' . 50 84-1 . 1239. 1 : « 5303o 0 -.4169.-1 .1644, 10 

a 04-1 .1472. 1. ;3125. 0 ,-.4202.-1 .1647. 10 , 
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